(N 

o 



in 



Relativistic nucleon optical potentials with isospin dependence in Dirac Brueckner 

Hartree-Fock approach 

Ruirui Xu and Zhongyu Ma 
China Institute of Atomic Energy, P.O. Box 275(41), Beijing 102413, China 

E. N. E. van Dalen and H. Muther 

Institut fiir Theoretische Physik, Universitat Tiibingen, 

Auf der Morgenstelle I4, D-72076 Tiibingen, Germany 

(Dated: February 27, 2013) 

The relativistic optical model potential (OMP) for nucleon-nucleus scattering is investigated in 
the framework of Dirac-Brueckner-Hartree-Fock (DBHF) approach using the Bonn-B One-Boson- 
Exchange potential for the bare nucleon-nucleon interaction. Both real and imaginary parts of 
isospin-dependent nucleon self-energies in nuclear medium are derived from the DBHF approach 
based on the projection techniques within the subtracted T-matrix representation. The Dirac po- 
tentials as well as the corresponding Schrodinger equivalent potentials are evaluated. An improved 
local density approximation is employed in this analysis, where a range parameter is included to 
account for a finite-range correction of the nucleon-nucleon interaction. As an example the total 
cross sections, differential elastic scattering cross sections, analyzing powers for n, p + '^'^Al at in- 
cident energy 100 keV ^ -E ^ 250 MeV are calculated. The results derived from this microscopic 
approach of the OMP are compared to the experimental data, as well as the results obtained with 
a phenomenological OMP. A good agreement between the theoretical results and the measurements 
can be achieved for all incident energies using a constant value for the range parameter. 
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I. INTRODUCTION 

The Optical Model potential (OMP) is one of the es- 
sential tools of exploring the physics of nuclear reac- 
tions. The investigations of OMP based on a fundamen- 
tal microscopic theory attracts a lot of interest because 
it provides a reliable basis to explore nuclear reaction 
mechanisms in particular for cases without experimental 
data. With the development of radioactive beam facili- 
ties, many new physical phenomena have been and will 
be observed for exotic nuclei and unstable isotopes be- 
yond the /3 stable line. A number of nuclei with large 
neutron excess will be produced and their structure and 
collisions are under investigation. Thus, the isospin de- 
pendence of the optical potentials becomes an important 
issue in the systematic study of OMP. However, this fea- 
ture was underdeveloped for lack of direct guidance from 
measurements in the past. Therefore, the main purpose 
of this work is to construct the nucleon OMP in a mi- 
croscopic way and discuss its features, especially with 
respect to the isospin dependence. 

Different models have been used to determine a mi- 
croscopic OMP [l|-13- They can be categorized as semi- 
phenomenological and ab initio approaches according to 
the nucleon-nucleon interaction and the approximation 
scheme used to describe the many-body systems. The 
semi-phcnomenological approaches normally begin with 
a phenomenological nucleon-nucleon (NN) effective in- 
teraction treated within a mean field approximation. 
Typical examples are the Skyrme-Hartree-Fock (SHF)[3| 
and the relativistic mean- field (RMF)0]. In these ap- 
proaches, the TV TV effective interactions are phenomeno- 



logically adjusted to reproduce nuclear matter saturation 
properties and the ground state properties of a few well 
known stable nuclei. The imaginary parts of OMP are 
determined in an independent manner through the eval- 
uation of the polarization diagram. On the other hand, 
the ab initio method, such as the Brueckner-Hartree- 
Fock (BHF)[l|,[ 3 , M ,[g| and Dirac-Brueckner-Hartree-Fock 
(DBHF) [a. [l0l - [l5| approaches, starts from a microscopic 
model of the bare NN interaction and the nuclear many- 
body problem is handled in a more sophisticated manner. 
Since the NN interaction is adjusted to describe the two- 
nucleon data, the many-body calculation docs not con- 
tain any adjustable parameter. Therefore such studies 
should have a larger predictive power when applied in 
the study of exotic nuclear systems. 

It is well known that the non-rclativistic BHF fails to 
reproduce the empirical saturation properties of nuclear 
matter without an extra three-body forcefJQ- An en- 
couraging success has been obtained when the problem is 
addressed in a relativistic framework, namely the DBHF 
approach[lO 12,16]. The saturation properties of nuclear 
matter in DBHF are described rather accurately without 
any need to add a three-body force. Moreover, another 
important feature of the relativistic approach applied to 
finite nuclei is, that the spin-orbit potential arises nat- 
urally from the coherent sum of the contribution from 
the scalar and vector potentials 0, which implies that 
this important ingredient of OMP occurs in the DBHF 
treatment without the need of any adjustments. 

The optical potential of a nucleon in the nuclear 
medium corresponds to the nucleon self-energy [17|. In 
the DBHF approach one can distinguish between two fre- 



qucntly used schemes [T^ to determine the relativistic 
components of the self-energy: the fit method versus the 
projection technique method. In the fit method, orig- 
inahy proposed by Brockmann and Machleidt [il|, the 
scalar and vector self-energy components are directly ex- 
tracted from the single particle energy. A fit to describe 
the momentum dependence of the single-particle energy 
in terms of constant scalar and vector components defines 
these self-energy components for each density considered. 
An attempt, which tries to extend this method in or- 
der to extract momentum dependent scalar and vector 
components by elaborate fitting procedures [l^l, suffers 
from large uncertainties. Therefore, only mean values 
for the self-energy components, where the explicit mo- 
mentum dependence has already been averaged out, can 
be obtained. In symmetric nuclear matter this method 
turned out to be rather reliable and predicts a relativis- 
tic decomposition of the self-energy in close agreement 
to the analysis of the projection technique. The exten- 
sion of this simple method to the isospin asymmetric 
nuclear matter, however, fails to determine the correct 
behavior of the isospin dependence of the nucleon self- 
energies [la, [2^. Therefore, this method is not very well 
suited for the construction of the microscopic OMP. 

In the other method, the projection technique method, 
the scalar and vector components of the self-energies are 
directly determined from the projection onto Lorentz 
invariant amplitudes. These projection techniques are 
rather involved but more accurate. It requires the knowl- 
edge of the Lorentz structure of the positive-energy- 
projected in-medium on-shell T-matrix. However, am- 
biguities [21[ arise due to the restriction to positive en- 
ergy states, since pseudo-scalar (ps) and pseudo-vector 
{pv) components can not uniquely be disentangled for on- 
shell scattering. As a consequence an unphysical strong 
momentum dependence of the nucleon self-energy is ob- 
tained, when using the ps representation scheme [22| . 
This strong momentum dependence mainly originates 
from the contribution of the single tt exchange. There- 
fore, it was suggested to separate the single meson ex- 
change contributions, i.e. the bare interaction V, and the 
high-order corrections AT [22] . The contributions of sin- 
gle TT and 1] exchange are in this approach treated with the 
complete pseudo- vector (pv) representation, whereas the 
pseudo-scalar (ps) representation is applied only on the 
remaining part of T-matrix, i.e. the high-order correc- 
tions and the single meson exchanges of the other mesons. 
This representation scheme for the T-matrix is called the 
subtracted T-matrix (STM). The application of this rep- 
resentation scheme to the single meson exchanges repro- 
duces the Hartree-Fock nucleon self-energies. Therefore, 
it is in princi ple equivalent to the method used in the 
work of Ref. [1^1, which not only considered symmet- 
ric nuclear matter but also pure neutron matter. This 
method was successfully applied to the investigation of 
the microscopic OMPJg as well as the nucleus-nucleus 
scattering with a double folding methodj23|. 

Recently an obvious improvement has been achieved 



in the DBHF calculations using projection techniques 
by the extension to the general case of isospin asym- 
metric nuclear matter, in which the neutrons and pro- 
tons are occupying different Fermi spheres jl3l - KLa |. As 
a consequence, one has different effective masses and 
self-energies for neutrons and protons. Furthermore, one 
has to deal with three different Pauli operators and in- 
medium interactions of the nuclcons, to be specific the 
nri, pp, and np ones. In contrast to the five independent 
helicity matrix elements for identical particles, six helic- 
ity matrix elements are independent in the np channel. 
Hence, instead of the five Lorentz invariants in the nn 
and pp channel, the T-matrix for the np channel is pro- 
jected onto six Lorentz invariants. With this method one 
can obtain the real part and the imaginary part of the 
self-energy for symmetric and asymmetric nuclear mat- 
ter. 

This method has been applied in this work to investi- 
gate the isospin-dependent relativistic microscopical op- 
tical potential (RMOP). To obtain the OMP for finite 
nuclei, an improved local density approximation (ILDA) 
method [l| is employed to determine the spacial distri- 
butions of the corresponding self-energies from the rele- 
vant density distribution and asymmetry of the finite nu- 
cleus under consideration. As an example we analyze the 
elastic scattering reactions of n,p-|-^^Al in terms of the 
RMOP and compare with experimental data and those 
obtained from an analysis with a phenomenological OMP. 

The paper is arranged as follows. The formalism of 
DBHF in STM representation is brieffy introduced and 
the nuclear matter properties with Bonn-B are described 
in Sec. II. In Sec. HI we present the isospin-dependent 
RMOP and some relevant discussions. The applications 
of RMOP to the neutron and proton elastic scattering off 
^^Al are shown in Sec. IV. Finally, a brief summary is 
presented in Sec. V. 



II. DBHF APPROACH IN ASYMMETRIC 
NUCLEAR MATTER 

A. DBHF approach 

In the relativistic DBHF scheme, the interaction be- 
tween nuclcons in nuclear matter is determined through 
the ladder approximation of the relativistic Bethe- 
Salpeter (BS) equation. 



T = V + i / VQGGT, 



(1) 



where T represents the nucleon-nucleon interaction ma- 
trix in the nuclear medium and V is the bare NN in- 
teraction, respectively. The Pauli exclusion principle is 
included by Q operator and the in-medium nucleon prop- 
agation is introduced by the Green's function G, which 
fulfills the Dyson equation. 



G — Go + GqTjG. 



(2) 



Go denotes the free nucleoli propagator, and the self- 
energy term E can easily be derived in first order through 
the following standard treatment 
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in which the self-energy contains the direct and exchange 
forms at the same time, and the momentum is integrated 
within the Fermi sea. Because the Eqs. (I!])-© are 
strongly coupled, they have to be solved iteratively until 
the convergence is reached. 

A one boson exchange potential (OBEP) model of the 
bare NN interaction is normally employed for the V term 
in Eq. ([T]), and Bonn-B NN interaction is selected in this 
discussion. Six mesons with different (J'^,T) are involved 
in the Bonn-B interaction, including two scalars a (0"'",0) 
and 5 (0+,l), two vectors w (1~,0) and p (1~,1), and two 
pseudo- vectors (pv) -q (0~,0) and tt (0~,1). For the de- 
scription of isospin asymmetric nuclear matter, one needs 
to distinguish between protons and neutrons. Therefore, 
this Bonn-potential code has been made suitable to treat 
distinct particles in the medium [ij]. As a consequence, 
one has three different in-medium interactions of the nu- 
cleons, to be specific the nn, pp, and np ones. The np 
channel has an additional independent amplitude com- 
pared to the five independent amplitudes in the nn and 
pp channel. 

The subtracted T-matrix (STM) representation has 
been applied for the projection of the T-matrix in the BS 
equation. The contributions of single tt and 77 are in this 
approach treated with the complete pv representation. 
The ps representation is applied only on the remaining 
part of the T-matrix, 



T,. 



T - V^^r, 



(4) 



In this STM representation the ambiguity problem of the 
self-energy from the contributions of tt and 77 mesons in 
isospin asymmetric nuclear matter can be minimized sat- 
isfactorily. Thorough discussions on generating the nu- 
cleon self-energies and comparing the reproduced bulk 
properties of the nuclear matter with the empirical val- 
ues for the symmetric and asymmetric matter can be 
found in Ref. [IMllll. 

The relativistic nucleon self-energy in the asymmet- 
ric nuclear matter is defined to satisfy the translational 
invariance of the homogeneous system, and its common 
Lorentz structure can be expressed as. 



E*(fc,fcF,/3) = E*(fc,fci.,/3) 



-7 



-7oEo(fc,fci^,/3) 



(5) 



where E^ is the scalar part of self-energy. Eg and E„ de- 
note the timelike and spacelike terms of the vector part, 
respectively. The superscript t is used to mark the proton 
and neutron, since one has to distinguish between neu- 
trons and protons in isospin asymmetric nuclear matter. 
It can be seen that the self-energies are the functions of 
the nucleon momentum (fc), density or Fermi momentum 
(fci?), and asymmetry parameter /? = (pn-pp)/ p, where 
Pn , Pp and p indicate the neutron, proton and total den- 
sities in nuclear matter. 



B. Nucleon self-energies and nuclear matter 
properties 

In this subsection, we present the bulk properties of 
nuclear matter as well as the nucleon self-energies in the 
DBHF approach with Bonn-B A^ A^ interaction calculated 
according to the procedure described in Ref. [l^]- The 
binding energies per nucleon (E/A) are shown in Fig. [1] 
as a function of the density, to be precise as a function 
of the corresponding Fermi momentum in symmetric nu- 
clear matter. Various values for the asymmetry param- 
eter f3 are considered in the range from 0.0 (symmetric 
matter) to 1.0 (neutron matter). The energy versus den- 
sity curves do not exhibit a minimum for /3 > 0.8. The 
saturation point calculated in symmetric nuclear matter 
occurs at kp = 1.33 fm'^ with E/A = -14.71 MeV. This 
is in good agreement with the empirical values kp = 1.36 
± 0.06 fm-i with E/A = -16 ± 1 MeV. Furthermore 
this figure demonstrates, that, at a given density, isospin 
asymmetric nuclear matter gets monotonically less at- 
tractive with increasing /3. 
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FIG. 1: Binding energy per nucleon in nuclear matter. 

For the example of symmetric nuclear matter at a 
Fermi momentum kp=1.35 fm~^ the real and imaginary 
parts of the nucleon self-energy components E^, Eo, and 
E^ are shown in Fig. [2] and Fig. [3l respectively, as func- 
tions of the on-shcll nucleon energy E. This single par- 



tide energy E is determined by 



^ = 



k^{l 



I]„)' + (M + E«)2-So-M. 



(6) 



In these figures, one can observe that the real parts of E^ 
and So change weakly with E in Fig. [2] especially for the 
energies above 100 MeV. A stronger energy dependence 
is obtained for the imaginary parts as can be seen in Fig. 
[31 Note that the imaginary parts vanish at the energy E 
which corresponds to the Fermi energy. The real space- 
like part fcEt, is rather small in comparison to the other 
two components, while the imaginary k'Ey turns out to 
be comparable in size with ImE^, ImEg. 

The equivalent potential to be used in a Schrodinger 
equation may in a first approximation (more detailed 
discussion see below) be represented by the difference 
Es — Eo- Hence the real parts of the self-energy compo- 
nents compensate each other to a large extent leading to 
a constant value for Eg — Eg of approximately -60 MeV for 
large energies E. This kind of cancellation between scalar 
and vector components is not observed for the imaginary 
part. Nevertheless also the imaginary part of E^ — Eq is 
essentially constant around -20 MeV at larger energies. 
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FIG. 2: Real parts of the nucleon self-energies at fcF=1.35 
fm"^ in symmetric nuclear matter. Note that k is the mo- 
mentum of single particle. 
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FIG. 3: Same as Fig[2] except for the imaginary parts. 



and £ = i? + M is the single particle energy, E is the 
kinetic energy of the nucleon in the free space and M 
indicates the mass of the nucleon. 

The Dirac potentials \J s and [/q with various values for 
the isospin asymmetries /? = 0.0, 0.2, 0.6, 1.0 at kp =1.35 
fm~^ are plotted as functions of energies in Fig. |4l Over- 
all, the strengths of both real and imaginary parts of 
scalar and vector potentials are sensitive to the asymme- 
try parameter in nuclear matter. The energy dependence 
is visibly infiuenced by the asymmetries for all Dirac po- 
tentials below energy E =100 MeV. The imaginary parts 
of the scalar and vector potentials are much weaker than 
those of the real parts. The behaviors of the imaginary 
parts are different from those of the real parts and their 
energy dependence is obviously stronger. The asymme- 
try dependence of the imaginary potentials for neutrons 
are much smaller than those for protons near the Fermi 
surface, but comparable at higher energies. 

The optical potentials of a nucleon scattering off finite 
nuclei are obtained by means of the local density approx- 
imation (LDA) as usual, in which the spacial function 
for RMOP at the incident energy E can directly be re- 
lated to the density (p), momentum (fc) and asymmetry 
parameter (/3) of nuclear matter by 



III. RELATIVISTIC MICROSCOPIC OPTICAL 
POTENTIAL IN FINITE NUCLEI 



C/LDA(r,S) = f/jVA/(fc,p(r),/3), 



(9) 



The Dirac equation of a nucleon in the nuclear medium 
can be written 



[d • p -f 7o(M+ [/*) + C/^] ** = £** , 
where ?7* and [/q are scalar and vector potentials, 
Ei-Ef,M ,,, -E*+eE* 
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(8) 



where Ulda and Unm correspond to the potentials of a 
finite nucleus and the nuclear matter, respectively, p{r) is 
the density at nuclear radius r and related to the Fermi 
momentum kp hy p = (2A;|,)/(37r^). To obtain the scat- 
tering amplitude and evaluate the observables for finite 
nuclei, a Schrodinger type equation is obtained by elim- 
inating the lower components of the Dirac spinor in a 
standard way. The equation for the upper components 
of the wave function is transformed into 
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FIG. 4: (Color online) Comparison of Us,Uo at ^^ = 1-35 fm^^. The real parts (upper panels) and imaginary parts (lower 
panels) are plotted for various /3 values. The black curves are the results for the symmetric nuclear matter, the blue curves 
denote those for proton and red curves for neutron. 
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where V^^t, T/*„. and F^„^„i„ represent the Schrodinger 
equivalent central, spin-orbit and Darwin potentials, re- 
spectively. These potentials are related to the scalar Us 
and vector Uq potentials by 
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where Vc is the Coulomb potential for a charged particle, 
D denotes a quantity defined as 



D\r) = A/ + e + Ul{r) - U^{r) - Vc. 



(12) 



One can see again that the spin-orbit term, Vg.o., emerges 
in a very natural way from the relativistic scheme. 

Based on the Dirac potentials above, the Schrodinger 
equivalent potentials can be computed following Eqs. 
dnil-lini). Firstly, in Fig. [H the central potential Vcent 
in the case of symmetric matter is shown as a function 
of incident energy at various densities, corresponding to 
kp^l-Ol, 1.14, 1.24, 1.33, 1.41 fm^^. For this figure we 
dropped the Coulomb contribution to Vcent- It is noted 
that the depths of real Vcent decrease with increasing en- 
ergy, and there is a crossing of curves for different densi- 
ties around E = 200 MeV. This implies that the so-called 
"wine-bottle bottom" shape appears at the the nuclear 
surface of finite nucleiQ. At energies below 50 MeV the 
imaginary potentials Vcent are deepest around kp = 1 
fm^^ and become weaker at higher densities. The shape 
of the imaginary potentials for finite nuclei contains a 
local maximum at the nuclear surface energies E' < 50 
MeV, which indicates a surface absorption. The situation 
is different at energies above 75 MeV, where the absolute 
value of the imaginary potential increases monotonically 
with the density, which exhibits the behavior of volume 
absorption. 

Additionally we compare the Vcent of proton and neu- 
tron versus the projectile energy E for various /? values 



at the fixed density kp ~ 1.35 fm~^ in Fig. [6l One finds 
that the depths of both real and imaginary potentials for 
neutrons are smaller than those for protons below £'=250 
MeV. We specially analyze the dependence of Vcent on 
the isospin asymmetry at several energies in Fig. [T] 
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FIG. 7: (Color online) The central potentials as functions of 
asymmetry /3 at kp — 1.35 fni~^. The black, red, green and 
blue ones correspond to the incident energies E = 10, 60, 
150, 250 MeV, respectively. Solid curves are for neutron and 
dashed ones for proton. 



FIG. 5: (color online)The Schrodinger equivalent central 
potentials as functions of the incident energies in symmet- 
ric nuclear matter at the nuclear Fermi momenta kp = 
1.01, 1.14, 1.24, 1.33, 1.41 fm~^. The real and imaginary parts 
of Vcent are plotted in left and right panels, respectively. 



lively taken as 0.037, 0.297 and 0.447. 
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FIG. 6: (color online) The Schrodinger equivalent central po- 
tentials as functions of the incident energies in the asymmetric 
nuclear matter at kp = 1.35 fm~^. The real part of Vcent is 
in the left panel. The right panel shows the imaginary part. 
The color notations are the same as Fig. [l] 

As an example for the applications of RMOP we study 
the experimental scattering data for the nucleon scat- 
tering off aluminum targets. An empirical formula by 
Negele [2J| is employed to generate the nucleon density 
distribution of the stable ^^ Al and two unstable neutron- 
rich isotopes •^'^'^'^Al, as shown in Fig. |51 The asymmetry 
is assumed constant in the interior of nucleus and deter- 
mined by /3 = (iV - Z)/A, where iV, Z, A indicate the 
numbers of neutron, proton and total nucleons of a finite 
nucleus. Thus, the asymmetries of ^7,37,47^^ ^^^ respcc- 




radii (fm) 
FIG. 8: Nucleon densities of ^'Al and two unstable isotopes 

37.47^1 

In this discussion, we adopt the finite range correction 
with Gaussian form to make modifications in the poten- 
tials via LDA approach, namely improved local density 
approximation (ILDA)[l[. The potential Uilda of the 
finite nucleus is finally obtained through the following 
integral, 

UiLDA{r,E) = {tV^)-^x (13) 

ULDA{r',E)expHr-7^'\yt^)d^r' , 



where t is the parameter, that represents the effective 
range of Ulda at radius r' . This is the only adjusted 
parameter in this calculation and is selected as a constant 
1.4 fm in all cases. 

We derive the Schrodinger equivalent potentials, Vcent ^ 



V* vh 



for n,p+^^'^^'^^Al through this ILDA ap- 



proach. The potentials of neutron scattered from the 
stable 27A1 at ^ = 10, 150, 250 MeV are plotted in Fig. 
m The depths of the real potentials decrease as the en- 
ergy increases and the imaginary parts change from the 
surface absorption to the volume absorption. These be- 
haviors are consistent with the phenomenological OMP. 
It can also be observed that the main contribution of 
RMOP comes from the central potential Vcent , the spin- 
orbit Vs.o. and Darwin Voarwin terms are the necessary 
supplements to the central potential. In Fig. [TUJ we con- 
struct the Vcent , for neutron scattered from ^^ Al and the 
neutron-rich '^^'^^Al. Due to the varying nucleon den- 
sity distributions and asymmetries, the depths of Vcent 
of the Al isotopes become smaller in the interior of the 
nuclei (r < 3 fm) with the increasing neutron numbers, 
while the potentials of heavier nuclei are expanded for a 
large radius. Meanwhile, we analyze Vcent of n+^^Al at 
E = 150 MeV in Fig. [TT] to check the impact of isospin 
asymmetry on the potential. It is found that the central 
potentials of ^^ Al are clearly affected by the isospin inside 
the nucleus. The impact is believed to become stronger 
for those nuclei with larger values for (3. 



IV. RESULTS AND DISCUSSION 

In this section, we discuss the angular distributions 
and analyzing powers of n,p-t-^^Al and the cross sections 
of 27Al(n,tot) at 100 keV< E <250 MeV. At very low 
energies the elastic scattering from the compound nu- 
cleus is sizable in addition to those from the OMP. In 
order to compare with the experimental data the Hauser- 
Feshbach model is utilized to determine the contribu- 
tion from the compound nucleus. For a particular in- 
cident energy, there are six open single particle emis- 
sion channels from the compound nucleus, including the 
neutron, proton, deuteron, tritium, a, Helium-S chan- 
nels. The first several discrete states are taken as the 
competing channels. For higher excitation energies, a 
continuum described by the GilbertCCameron level den- 
sity formula|25| is provided as compensation. Therefore, 
the computed observables contain the contributions from 
both shape and compound elastic scattering processes. 
All calculations mentioned above are performed with the 
assistance of the APMN code [26}. 

The present calculations are systematically compared 
with the experimental data and the results obtained with 
the global phenomenological OMP, Koning-Delaroche 
(KD) potential ^ in Figs. [T^lfTSl The experimental 
data in these figures, the total cross sections CTfot, the 
differential elastic scattering cross sections da/dil and 
the analyzing powers Ay (9), are assembled in Table U- 
im The calculated observables of the neutron and pro- 
ton induced reactions are both in good agreement with 
experimental data. Moreover, the results of the RMOP 
are very similar to those with the global KD potential. 
In some energy regions our predictions are even better 
than those with the phenomenological ones, such as the 



TABLE I: atot, da/dQ and Ay{9) database for n-l-^'^Al 



Data 



Rcf. 



Encrgy(MeV) 



O'tot 



da/dO, 



Ay{9) 



[28J 


0.25 - 27.0 


[29] 


5.3 - 250.0 


[30] 


2.0 - 80.1 


m 


3.2 


[32] 


5.4, 6.4, 7.5, 8.6 


[33] 


7.6 


m 


10.9, 13.9, 16.9 


[35] 


15.4 


[36] 


18.0, 20.0, 22.0, 25.0, 26.0 


[37] 


84.0 


[38] 


96.0 


[39] 


136.0 


[33] 


7.6 


[40] 


14.0, 17.0 


[35] 


15.4 



TABLE II: da/dO. and Ay{e) database for p-|-^^AI 



Data 


Rcf. 


Energy(MeV) 


da/dO. 


4^ 


0.8, 1.0, 2.0, 3.0 




42] 


17.0 




43] 


17.5 




4^ 


28.0 




45] 


61.4 




46] 


92.9, 95.7 




43 


142.0 




48] 


156.0 




[49, OT] 


160.0, 177.0, 183.0 




[5^ 


185.0 


Ay{e) 


[52] 


13.7, 15.0, 17.0, 19.0, 21.0, 23.0, 25.0 



angular distributions of p-|-^^Al a.i E > 95.7 MeV in Fig. 
[T3l As a conclusion, it is proved that the present RMOP 
based on the DBHF approach is suitable to describe the 
potentials of n,p-|-^^Al, especially for the spin-orbit term 
Vs.o., which is validated through the apparent agreement 
to the measured analyzing power data in Fig. 1141 

We particularly investigate the DBHF central potential 
of n+2'^Al at ^ = 150 MeV with the ILDA approach and 
the conventional LDA approach in Fig. [161 The Vcent 
of the phenomenological KD model is also displayed in 
this figure. It is observed that the depths of the real and 
imaginary Vcent in LDA method are obviously deeper and 
exhibit a smaller range than those of the KD potentials. 
This differences are reduced after the finite range correc- 
tions in ILDA approach are taken into account, and the 
corresponding RMOP looks closer to KD ones. There- 
fore, rather similar scattering results in Figs. 11211151 are 
obtained with both potentials. 



V. SUMMARY 

In this work, we study the RMOP for finite nuclei in 
the framework of the DBHF approach with projection 




Radii (fm) 



Radii (fm) 



FIG. 9: (Color online) The Schrodinger equivalent central, spin-orbit and Darwin potentials as functions of radius for n+^^AI 
at E =10, 150, 250 MeV. 
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FIG. 10: The Vcent of n-f^^'^^'-'^Al at £ = 150 MeV. The 
real parts are shown in the left and the imaginary ones in the 
right. 
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FIG. 11: The Vcent of n-l-^^Al in the isospin dependent mode 
and those for neglecting the isospin asymmetry at _E = 150 
MeV. The left panel shows the real parts, the right one for 
the imaginary parts. 



techniques using the Bonn-B bare NN interactions. Spe- 
cial attention is paid to the isospin dependence. The 
subtracted T-matrix representation is applied to mini- 
mize the ambiguities, which arise due to the restriction to 
positive energy states in the DBHF approach, in the de- 
termination of the nucleon self-energies for isospin asym- 
metric nuclear matter. In this way the precise form of 
the nucleon self-energies with its momentum and isospin 
dependence are determined. Thus, the nucleon RMOP 
is obtained from the results of these DBHF calculations. 
An ILDA method is adopted to intimately relate the den- 
sity and isospin dependence of the RMOP in nuclear mat- 



ter to the radial dependence in finite nuclei. A widely 
used empirical density distribution is employed for finite 
nuclei. Only one free parameter in ILDA is adapted in 
RMOP, which is fixed to be t = 1.4 fm for all calculations. 
As an example of the applications we investigate the 
elastic scattering of n.-p^'^'^ K\ at 100 keV < E < 250 
MeV. The predicted observables in this approach are 
compared with experimental data and with results from 
phcnomcnological KD potentials. Good agreement is ob- 
tained. In addition, we perform a first analysis of the 
RMOP for the unstable neutron-rich nuclei '^^^^^Al. It is 
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FIG. 12; Angular distributions ol n+'^^Al elastic scattering 
at 3.2 MeV < £ < 136 MeV. The curves and data points at 
the bottom represent their true values, the others are multi- 
plied by factors of 10 , 10 , etc. The solid curves and dashed 
ones indicate the results calculated with RMOP and KD po- 
tentials, respectively. The dots are experimental data hsted 
in Table U 
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FIG. 13: Same as in Fig |12l except for p+'^'^Al elastic scatter- 
ing for induced energy 800keV< E < 185.0 MeV. The experi- 
mental data are shown in Table HIl 



demonstrated that the isospin asymmetry has an obvious 
impact on the OMP potentials. 

To apply this RMOP in a convenient way and make it 
accessible to other research groups, explicit expressions 
for the nucleon scalar and vector potentials might be de- 
sired. In this work, we systematically analyze the en- 
ergy, density and isospin dependence of the nucleon self- 
energies and the Dirac potentials in nuclear matter and 
the Schrodinger equivalent potentials in finite nuclei. It 
is observed that the potential is sensitive to the isospin 
parameter /S, but some features related to the energy 



dependence arc not sensitive to /?. The full parameteri- 
zation of this RMOP and the extensive application will 

be studied in the following work. 
15 
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FIG. 14: Analyzing powers of p,n-|-'^^Al elastic scattering. 
The curves and data points at the bottom represent their 
true values, the others are shifted by increments of 2, 4, etc. 
The notations are the same as in Fig. 1121 The experimental 
data are listed in Table IHIll 
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E (MeV) 

FIG. 15: The cross sections of n-|-'^'^AI(n,tot). The experi- 
mental data are listed in TableU] The Phenomenological OMP 
indicates the results calculated with the KD potentials. 
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